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We make a matrix-product variational approach to spin-i ferromagnetic-antiferromagnetic bond- 
alternating chains with anisotropy on their ferromagnetic bonds, especially under the open boundary 
condition. The rich phase diagram containing the Haldane, large-D, and two types of Neel phases 
is well reproduced with only two variational parameters. The on-bond anisotropy has a significant 
effect on the ferromagnetic coupling between neighboring spins and induces novel edge states peculiar 
to spin-i chains. 

PACS numbers: 75.10.Jm, 75.40.Cx, 02.70.Wz, 75.40.Mg 



Haldane's conjecture 0, Q sparked renewed interest 
in low-dimensional quantum magnets and led to ex- 
tensive explorations of spin gaps— energy gaps in mag- 
netic excitation spectra. While the early research cen- 
tered on spin-1 antiferromagnetic chains, spin-i bond- 
alternating chains have been attracting further inter- 
est due to their enriched ground-state properties 0, 0, 
0, 0, 0, I3- A spin-i Heisenberg chain with alternat- 
ing ferromagnetic and antiferromagnetic couplings con- 
verges to the spin-1 antiferromagnetic Heisenberg chain 
as the ferromagnetic coupling tends to infinity, where 
the ground state remains unique and gapful for all the 
ratios of the two exchange couplings [2|. Thus the 
spin-1 Haldane gap turns out to be continuously con- 
nected to the gap originating in decoupled singlet dimcrs. 
Such a consideration can be verified for existent bond- 
alternating chain compounds such as IPACuCls (IPA = 
isopropylammonium = (CH3)2CHNH3) and (4- 
BzpipdH)CuCl3 (4-BzpipdH = 4-benzylpiperidinium = 
C12H18N) (Iflj. These materials indeed reproduce many 
of observations common to Haldane-gap antiferromag- 
nets With increasing temperature, the 

,tures disapi 
magnetic behavior [isl Il7l | . 

Apart from paramagnetic spin ^ ^, ^^^^^^^.^^ ^j,,^^^ ^ 
degrees of freedom lie in spin-1 Haldane-gap antiferro- 
magnets El E El, 22, 23, ^ 25, 26, 27, 2^. 
In more general, there appear fractional spin--^ degrees 
of freedom on boundaries of antiferromagnetic Heisen- 
ber g ch ains with integral spin S at low temperatures 
[29I l3Cl|| . The chain-end fractional spins are understand- 
able in view of the valence-bond-solid states l3lL l32i and 
were actually observed for both spin-1 [s^ l34l 135] and 
spin-2 [2^ antiferromagnetic chain compounds. Such ob- 
servations are peculiar to the Haldane phase. Hence both 
quantum and classical spin-i degrees of freedom may 
be observed for Haldane-gap materials composed of spin 
i's. Thus motivated, we study spin-i ferromagnetic- 
antiferromagnetic bond-alternating chains with particu- 
lar emphasis on their edge states, making good use of 
the matrix-product representation as well as a quantum 



effective spin-1 features disappear into the spin-i para- 



^'s, quantum spin-^ 



Monte Carlo method. 

The Hamiltonian of our interest is given by 
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where Jp and Jaf are both set positive. The on-bond 
anisotropy D, which is related to a single ion anisotropy 
in the case of spin 1, originates from possible dipole- 
dipole and/or anisotropic exchange interactions between 
the ferromagnetically coupled spin i's When we 

consider a wave function of the matrix-product type: 

\^) ^Tiigi(E)g2(E)---(E)gN], (2) 

employing the complete set of states on the jth. ferromag- 
netic bond: 
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the decoupled dimers for Jp ^ are described by 
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while the optimum ground state for Jp — s- 00 is given by 
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The g matrices (dJ, with their indices contracted, create 
singlet bonds in betwee n l37l and end up as the spin-1 
valence-bond-solid state [32|- The g matrices Q^can be 
used as the basis for a variational calculation [38| of the 
Hamiltonian by allowing different amplitudes for the 
singlet and triplet contributions and reducing the rota- 
tional invariance as 
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TABLE L The matrix-product variational (MP) and quantum 
Monte Carlo (QMC) calculations of the ground-state energy 
per site for the isotropic {D — 0) chain. 



Jf/Ja 



MP 



QMC 



0.2 
0.5 
1.0 
2.0 
5.0 



-0.37678 
-0.38538 
-0.41197 
-0.49435 
-0.82583 



-0.3768(1) 
-0.3855(1) 
-0.4125(1) 
-0.4976(1) 
-0.8426(1) 



The isotropic chain is described by equaUzing c to a and 
the two extreme cases Q and (jSJ are indeed included 
in this ansatz. Then, eq. ^ gives a variational ground 
state of the periodic chain, while each of the four elements 
of the 2x2 matrix gi ^ g2 ^ ■ ■ ■ ^ 9n corresponds to a 
ground state of the open chain with fixed spins at each 
end. In the Haldane phase, the four states of the open 
chain are quasi-degenerate 39]. 

The variational energy (*|7Y|*)/(vI'|vI/) 
culated as 
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and its bounds are compared with quantum Monte Carlo 
findings in Table HI The variational estimates are less 
precise in the region of strong ferromagnetic coupling 
but still agree with the numerical values within two- 
percent error even at Jf/Jaf = 5.0. The spin correlation 
function and local moments induced on a chain end are 
also calculated and verified in Fig. ^ Considering that 
the correlation length is considerably underestimated by 
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FIG. 1: The matrix-product variational (MP) and quantum 
Monte Carlo (QMC) calculations of the spin correlations (the 
upper three) and moments (the lower three) in z direction for 
the isotropic {D — 0) open chains with ferromagnetic cou- 
plings at each end. 
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FIG. 2: Quantum Monte Carlo calculations of the Curie com- 
ponent in the magnetic susceptibility as a function of temper- 
ature for the isotropic (D = 0) open chains with ferromagnetic 
(o) or antiferromagnetic (x) couplings at each end. The solid 
and dotted lines denote the Curie susceptibilities due to two 
spins of S = i and a single spin of 5 = 1, respectively. 



wave functions of the matrix-product type in general [40j , 
our variational findings are fairly good. The variational 
expression for the z-component correlation length ^, 



e 2 



|6p + |cp + 2|a|' 
\b\^ + \c\^-2H 
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is useful especially in the region of weak ferromagnetic 
coupling. The edge moments, which are assembled 
into an effective spin i, are characteristic of the Hal- 
dane phase 0, ll^l and should therefore be observed 
for the present system with nonmagnetic impurities as 
well. Besides direct observations, we can detect them 
through magnetic susceptibility measurements [l^ l3(Tj |. 
FigureElshows the Curie component of the susceptibility, 
X°2N+2 - X2n' = XC, where x^"' ^nd xT' ^^'^ suscep- 
tibilities of the periodic and open chains with L spins, 
respectively. Two excess moments of 5 = ^ are found at 
both low and high temperatures. The low-temperature 
ones are quantum mechanically correlated effective spins 
induced at each end, whereas the high-temperature ones 
are paramagnetic spins of classical aspect. In the region 
of strong ferromagnetic coupling, such a difference can be 
understood better, because the distinct features of para- 
magnetic spin I's appear at intermediate temperatures, 
where the Haldane phase is thermally broken but the fer- 
romagnetic exchange interactions still survive. The open 
chain with antiferromagnetic couplings at each end ex- 
hibits no edge moment at low temperatures even in the 
Haldane phase. 

Our variational scheme is still useful for the anistropic 
chain and well reproduces its rich phase diagram, as is 
shown in Fig. |3| With increasing \D\, the Haldane gap 
decreases to zero and there alternatively appear three 
distinct phases according to the sign of D and the ratio of 
Jf to Jaf- The variational wave function illuminates in 
itself how the spin correlations vary with the anisotropy. 
Moving across the phase transitions, we plot in Fig. 0] 
the optimum variational parameters. Considering that 
eq. © is rewritten as 
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FIG. 3: The matrix-product variational (MP) and numerical 
diagonalization (Exact) 4] calculations of the ground-state 
phase diagram. Schematical representations of each phase are 
presented for reference, where the arrow (the bullet symbol) 
and the segment denote a spin i with its fixed (unfixed) pro- 
jection value and a singlet pair, respectively, while the circle 
means an operation of constructing a spin 1 by symmetrizing 
the two spin i's inside. 



ment remains unchanged and monotonically grows into 
a long-range order, which is convincing on the analogy 
of our experience for a spin-1 chain |4l| . On the ways 
to the Neel I and large-Z? phases, on the other hand, 
we have novel observations peculiar to the spin-^ chain. 
Figure [SI shows the edge moments as functions of D. On 
both the ways from the Haldane to Neel I and large- 13 
phases, the spin alignment changes with increasing D 
from Ttiittii • • • to nnnn • • • via tttttttt . . .. The 
ferromagnetic alignment at a chain end is another inter- 
esting effect of D and is induced near D/Jaf = 0.22 
and D/Jaf = 0.5 in the cases of Jf/Jaf — 0.5 and 
Jf/Jaf — 2.0, respectively. The magnetic moments in 
z direction form a long-range order in the Neel phases, 
whereas they vanish in the large- 13 phase. The total mo- 



ment induced at a chain end, J2j=ii i^ij- 
is variationally expressed as 
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FIG. 4: The optimum variational parameters as functions of 
the ferromagnetic coupling and the on-bond anisotropy. 



where {\t-), ± i\ty),)/V2 = K±)„ \t^)j - |io)j, and 
CT° and (cr^, tJ*', cr^) are the 2x2 unit matrix and the 
Pauli matrices, respectively, we find that the anisotropy- 
induced phases are characterized in terms of S{S + 1) = 



iS2 



and S*^ 



On 



5|j- as i) Large-D: 



S' = 1, S*^ = 0; ii) Neel I : S*^ = 0; iii) Neel II: S = I, 
= ±1. Thus, the large-D and Neel II phases appear 
in a spin-1 chain as well, while the Neel I phase is charac- 
teristic of the spin-i chain. No ferromagnetic correlation 
between any neighboring spins in the Neel I phase can be 
understood by rewriting the Hamiltonian as 



, N 



2D — ,/f)'S'|,_iS'|,- + JAFS2j ■ S2j+l 



(10) 



With increasing D, the ferromagnetic coupling is reduced 
and turned into an antiferromagnetic one. 

We have already observed in Fig. ^ that the spin-i 
open chain also exhibits effective moments of iS = ^ at 
each end. In the isotropic case, there appears an align- 
ment of the type TTiiTTii ■ ■ ■ at a chain end. On the 
way from the Haldane to Neel II phases, such an align- 



S: 



edge 



/|5|2 + |c|2 -2|a|2x^/2 
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The absolute value of the fraction + |cp — 

2|ap)/(|6p + |cp + 2|ap) is smaller than unity in the Hal- 
dane phase, while it is fixed to unity in any other phase. 
Thus, the effective moment of 5 = ^ persists in the edge 
of an infinite chain throughout the Haldane phase and 
disappears with the collapse of the Haldane gap. 

Possible dipole-dipole and/or anisotropic exchange 
interactions between ferromagnetically coupled two 
spins act as a fictitious single-ion anisotropy jl5j | and 
high-frequency electron-spin-resonance experiments on 
IPACuCla J£j indeed revealed the nondegenerate triplet 
excitation. Therefore, we take more and more inter- 
est in nonmagnetic-ion substitution at the Cu site of 
spin-i ferromagnetic-antiferromagnetic bond-alternating 
chain compounds. The present system may be com- 
pared with further polymerized bond-alternating cop- 
per(II) complexes such as the trimerized chain com- 
pound 3CuCl2-2dx (dx = 1,4-dioxane = C4H8O2) El 
and the tetramerized chain compound Cu(3-Clpy)2(N.-^)^ 
(3-Clpy = 3-chloropyridine = C5CIH4N) [pL which be- 
have as a spin- 1 critical antiferromaenet |44l l45l| and a 
spin-(|, i) ferrimagnet 0|, respectively. With bond al- 
ternation, impurity-induced magnetic effects vary with 
the location of the dopant impurities in general, as is 
shown in Fig. El A comparative study on the edge states 
of various bond-alternating chain compounds will lead to 
brand-new observations. 

The authors are grateful to K. Hida, K. Okamoto and 
T. Sakai for valuable comments. This work was sup- 
ported by the Ministry of Education, Culture, Sports, 
Science and Technology of Japan, and the Iketani Sci- 
ence and Technology Foundation. 
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FIG. 5: The matrix-product variational (MP) calculations of the spin moments in z direction for the anisotropic open chains. 
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